clear
clc
x=0:0.1:6;
[x,y,z]=cylinder(9-(5/6)*x,50);
surf(x,y,z);
s1=16*pi;s2=81*pi;h=6;
V=1/3*(s1+s2+sqrt(s1*s2))*h